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The quark-gluon plasma formed in heavy ion collisions contains charged 
_ chiral fermions evolving in an external magnetic field. At finite density 

of electric charge or baryon number (resulting either from nuclear stop- 
ping or from fluctuations), the triangle anomaly induces in the plasma 
the Chiral Magnetic Wave (CMW). The CMW first induces a separation 

■ of the right and left chiral charges along the magnetic field; the result- 
S_i ' ing dipolar axial charge density in turn induces the oppositely directed 

vector charge currents leading to an electric quadrupole moment of the 
quark-gluon plasma. Boosted by the strong collective flow, the electric 
quadrupole moment translates into the charge dependence of the ellip- 
tic flow coefficients, so that V2(iv + ) < V2(tt~) (at positive net charge). 
Using the latest quantitative simulations of the produced magnetic field 
and solving the CMW equation, we make further quantitative estimates 
of the produced V2 splitting and its centrality dependence. We compare 
the results with the available experimental data. 
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1 Introduction 



In (off-central) heavy ion collisions, the moving ions produce a huge mag- 
netic field that could reveal an interesting dynamics in the plasma. The chi- 
ral anomaly present in quantum electrodynamics (QED) induces interesting 
"macroscopic" effects in the plasma: the magnetic field together with the pres- 
ence of a non- vanishing electric or axial charge density lead to axial or electric 
currents. These currents lead respectively to the so-called Chiral Separation 
Effect (CSE) [TJ[2J[3] and the Chiral Magnetic Effect (CME) H [5j 1 . 

The CME needs an initial axial charge fluctuation. It can be created by 
sphaleron transitions that are essentially random, i.e. the sign of the effect 
varies from event to event. Note that the QCD sphaleron rate is expected 
to be large (r - 0.3 T 4 ) and numerous sphaleron transitions happen in a 
single collision. The presence of strong magnetic field was found to increase the 
sphaleron rate [10], but for the realistic values of magnetic field in heavy ion 
collisions the change in the rate is small. If we suppose that sphalcrons create a 
non-zero axial chemical potential (ia, the QED anomaly induces a vector current 
jy = iprfip in the presence of an external magnetic field (CME): 

The vector current^ induces a movement of electric charges leading to an 
electic dipole in the plasma i.e. to "local parity violation" . Recently, STAR JT2l 
[T3"] . PHENIX P3] and ALICE [TS] collaborations observed charge asymmetry 
fluctuations possibly providing an evidence for the CME. The interpretation of 
the results however is still under intense discussion, see e.g. [TBI fTTl [T8] , 

The CSE requires an initial density of vector charges e.g. a finite baryon or 
electric charge number density, which can be present in the plasma not just as 
a fluctuation but as a non-zero background. Together with a magnetic field, the 
QED anomaly leads to an axial current: 

3 a = -r-^-nvB, (2) 

where /iy is the vector chemical potential. 

The coupling of the vector and axial currents induced by the anomaly in 
the presence of magnetic field brings to the existence a novel type of collective 
gaplcss excitation in plasma - the Chiral Magnetic Wave (CMW) [19]. While 
the CMW can exist even at zero density of vector charge, it leads to a very 
interesting phenomenon when the finite vector charge density is present [20] . 
One can visualize it in the following way: first the CSE leads to a dipolar 
separation of axial charges along the B field direction; then the finite density of 
axial charge in turn induces the oppositely directed CME currents resulting in a 
charge quadrupole that has the same sign in all collisions (as long as the density 

lr The existence of this current in the quark-gluon plasma has been established in lattice 
QCDxQED [□]. 
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Figure 1: Charged pair azimuthal correlations measured by the STAR at RHIC 
and the ALICE at LHC from Ref. [TS]. 



of vector charge is positive). Consequently there are more positive charges at 
the poles of almond-shape fireball (since B is primarily out-of-plane) than at 
the equator (in the reaction plane). This eventually gives rise to the difference 
in elliptic flows between positive charged particles and the negatively charged 
ones. In other words, the combination of CME and CSE, that is the chiral 
magnetic wave, leads to a charge dependence of the elliptic flow that survives 
even after averaging over events. 



2 The charge dependence of elliptic flow 

We will estimate here the size of the charge quadrupole created by the chiral 
magnetic wave in heavy ion collisions. 



2.1 Chiral magnetic wave 

The co-evolution of the vector and axial currents can be described in the frame- 
work of Chiral Magnetic Wave (CMW) equation [19]. Let us first give a short 
review of its derivation. We can rewrite the anomaly formulas (|1|2|) as 



j v \ _ N c eB ( 1 \ / n v 



3a) 2tt 2 \ 1 J \ [i a 
On the other hand, the chemical potential depends on the currents: 
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Remembering that the chemical potentials are derivatives of the free energy 

Wi 

we can define the susceptibility matrices a: 

Under parity transformation V — > —V and A — > A, so that ay a = &av = 0. 
Moreover one can show that ayy ~ &AA = « in the large N c limit. Inserting 
these results back into equation (j3J), we get 



(6) 



(7) 



j v \ N c eBa ( 1 \( j° 

3 A J 2tT 2 V 1 J V # 

We can rewrite these equations in the basis of left and right spinors: 

f N c eBa\ n , N 

3^ = ^{-^-)0U + - ( 8 ) 

The final step is to combine the equation ((HJ with the conservation laws d^j^ R = 
that are approximately valid in time scales much shorter than the quark 
mass; this is easily satisfied in real experiments. On general grounds one should 
also include diffusion, which arises at the next order in the gradient expansion. 
Denoting Xi the direction of the magnetic field, the evolution of the CMW is 
described by the wave equation 

(do T vdt - D L d\ - D T d 2 T ) f L>R = 0, (9) 

with v = the velocity of the wave and (Dt) the longitudinal (trans- 

verse) diffusion constant. Note that the constants a, Dl, Dt are B and T 
dependent and such that v < c. We use in the following the values of and 
D L that were computed in rcf. |19j in the framework of the Sakai-Sugimoto 
model in the large N c quenched approximation. 

To describe the initial "almond" configuration of the QCD matter we use 
the KLN model |H] based on parton saturation and kx factorization. Au — Au 
collisions have been simulated, with realistic Woods-Saxon nuclear densities. 
The axial chemical potentials at initial time are set to zero. 

We consider the propagation of the u and d quark currents in terms of 
two different CMWs. We omit here the strange quark as it has no net den- 
sity in the plasma. As the quark flavours have different charges g/, we use 
v(qfB), DL,r{qfB) as coefficients for the CMW. The total electric charge dis- 
tribution is obtained by summing the different waves as 

f: ^n(f;: J ■/;/)■ do) 
/ 

The resulting distribution is shown in fig. 12.11 We see that the excess of positive 
charges move towards the pole of the fireball whereas the center contains less 
positive charges, leading to a charge quadrupolc. 
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Figure 2: Initial baryon density and evolution of the excess of charges in the 
plane transverse to the beam axis. For clarity, we have subtracted the charge 
density distribution without the CMW in the second picture. Magnetic field 
strength used eB = 2.8m^., lifetime of magnetic field r = 4 fm, temperature 
T = 165 McV, impact parameter 6 = 3 fm. 



2.2 Prom charge quadrupole to elliptic flow 

In principle the CMW evolves in an expanding background. A full computation 
would require rewriting the CMW in an expanding medium and performing a 
full hydrodynamical simulation. Here for a first estimate of the effect, we will 
make the following approximations. We first compute the charge distribution 
of the wave evolution in a fixed background, and then include the effect of the 
flow. This would give precise results if the magnetic field has a short lifetime. 
In the case of a long lifetime, the propatation of the wave and the expansion of 
the medium will get entangled, which is a more complicated case that we will 
further study elsewhere. The dominant flow effect is the strong radial flow that 
correlates the emitted particle's transverse momentum magnitude and direction 
with its spatial radial position and azimuthal angle. We thus assume that the 
charge asymmetry distribution (at the end of CMW evolution) is simply carried 
over by the radial flow. 

After the evolution of the CMW, the azimuthal distribution of net charge 
density clearly bears quadrupole component, i.e. the quark-gluon plasma be- 
comes an electric quadrupole, sec Fig. 12.21 As mentioned above, the strong radial 
flow aligns the particles' momenta along the direction of the flow. Therefore the 
net charge detected at an angle <\> is proportional to 



The azimuthal anisotropy of net charge distribution is dominated by the elliptic 
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Figure 3: Schematic demonstration of charge quadrupole boosted by strong 
collective (radial) flow. 



component (i.e. the 2nd harmonics), and we approximately have 

d{N+ ^ N ~ ] = (N+ -N-)[l- r e cos(20)], (11) 

with N± the total number of ± charged particles and r e the ratio between the 
2nd harmonic of the net charge distribution to the total net charge: 

r e = 2^, (12) 



q e = - I rdrd4> cos(20 (j e °(r, 4>) ~ f e ,B=o(r, 4>)) 
p e = / rdrd4jl{r,4>). 



Note that in the definition of q e we focus on the charge quadrupole entirely due 
to the CMW evolution by subtracting out the net charge distribution without 
CMW. Supposing that the hydrodynamical expansion of the plasma is not per- 
turbed by the net charge distribution (which indeed is a very small perturbation 
to the overall bulk evolution), the total (charge independent) particle distribu- 
tion will still follow the usual elliptic flow pattern [1+2ij 2 cos(2</>)]. Using Eq. tjTTj) 
we can then write down the positive/negative charge distributions that include 
both the usual elliptic flow effect and the charge quadrupole contribution (up 
to 0(v 2 x r e A)): 

^±=N ± [1 + (2v 2 T r e A) cos(2<£)] , 

where A = is the net charge asymmetry. It is then obvious to see that 

the elliptic flow is charge-dependent vf = v 2 T r e A/2 with the splitting being 

V 2 ~ V t ~ r e A ■ (13) 

The splitting is linear in the net charge asymmetry A, with the slope r e deter- 
mined from the net charge distribution due to the CMW evolution. 
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2.3 Evolution of the charges after the plasma phase 

Upon hadronization, the electric charges carried by the quarks are transferred 
into pions, kaons, protons and other hadrons. However, various hadronic re- 
actions in the "afterburner" phase of the heavy ion collision may considerably 
alter the charge asymmetry pattern from the plasma phase. Protons and an- 
tiprotons as well as kaons and antikaons have quite different cross sections in the 
hadronic matter at finite baryon density, and this may well mask the effect one 
would like to look for. Positive and negative pions seem to be the best channel, 
as they have a relatively small difference in the absorption cross sections and 
are produced in large quantities. Therefore we suggest to measure the electric 
quadrupolc moment of the plasma by using the difference of the elliptic flows of 
pions, V2(tt~) — V2(tt + ) as the observable. 

3 Results and comparison with experiments 

As already discussed, the creation of the charge quadrupole from CMW requires 
an initial vector charge density. In relatively low-energy heavy ion collisions the 
plasma naturally acquires sizable initial baryonic charge density (and a relatively 
large charge asymmetry A) due to nuclear stopping effects. Therefore it is 
natural to first look for the predicted elliptic flow difference ^(tt ) — V2(tt + ) = 
r e A in low energy collisions, as proposed and estimated in [5D]. At high energy, 
A becomes small on average as the initial baryon density is small. However, 
as a result of event-by-event fluctuations, the charge asymmetry A of hadrons 
measured in a limited kincmatical domain (e.g. in a slice of rapidity) always 
fluctuates and can be quite sizable in a fraction of events. Therefore a second 
way to look for the predicted v 2 splitting is to measure the ^(tt - ) — v 2 {^ + ) 
versus the charge asymmetry A in high energy collision events (by binning them 
according to the charge asymmetry). In that way one can (1) test the predicted 
linear relation between the splitting and the charge asymmetry, and (2) extract 
the slope that can be directly compared with the quadrupole quantity r e from 
CMW computation. This very useful strategy has been recently developed and 
used by the STAR collaboration to measure the charge quadrupole at high 
collision energies [22] . 

3.1 Results from the numerical simulation 

In the original estimate |20| of the quadrupole moment, we assumed that the 
initial up and down contents of the plasma followed the up and down quark 
content of the gold ion. This was motivated by the picture of baryon stopping 
appropriate at low energies. Here we present slightly different results, where 
we have fixed the initial charge density and realize it with equal amount of up 
and down quarks. This is motivated by comparison with the STAR results at 
high energies where the amount of charge was measured directly. The resulting 
difference is however very small numerically and is well below other sources of 
incertitude. The main difference with the simulation of Rcf. [2D] is that we have 
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Figure 4: Magnetic field from recent event-by-event computations in Ref. [23) 
(BS) and Ref. [24] (DH). 



now used the state-of-art computation (on event- by- event basis) of the magnetic 
field (shown in Fig J3.1|) from Bzdak and Skokov (BS) [23], as well as from Deng 
and Huang (DH) [24] instead of the analytical estimate of magnetic field from 
Ref. [6] . Note that due to differences in their implementation, the results of BS 
are larger by a factor 2 than the results of DH, but both are significantly larger 
than the estimate of Ref. [5] for large impact parameters. 

The time evolution of the magnetic field in the plasma critically depends on 
the medium's transport properties. If the electric conductivity is large, mag- 
nctohydrodynamics is applicable and the magnetic flux is conserved. In the 
opposite limit, the plasma effects are negligible and the magnetic field decreases 
quickly as the spectator ions fly away. An estimate of the effect of electric 
conductivity on the lifetime of the plasma has been presented in |28| ; a more 
detailed study of this issue is clearly needed. Since the overall effect generated 
by the magnetic field would be integrated over time, instead of convoluting with 
the complicated time-dependent B field, we will compute the magnitude of the 
charge quadrupole by assuming the magnetic field to stay constant from to = 
till the end of the B field lifetime. The resulting electric quadrupole moment 
r e is shown in Fig. 13. li as function of the impact parameter b. 

3.2 Experimental results 

The first experimental evidences for the (global) electric quadrupole effect were 
found for 0-80% centrality events by the STAR collaboration at RHIC [25] and 
the CERES collaboration at SPS [26]. The sign and the magnitude of the 
difference in the positive/negative pion vi appear to agree with the prediction 
[20] . For instance for %fs = IIGeV collisions at RHIC, with A ~ 0.3 and 
r e ~ 0.01 for t = 4 fm using the magnetic field of (DH), one expects the 
difference in elliptic flow ~ r e A ~ 0.003 (with of 7r~ bigger than ir + ), which 
is 10% of the overall elliptic flow at this energy. These observations provide an 
important evidence for the proposed CMW effect. Note however that at such 
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Figure 6: The charge dependence of the elliptic flow from minimum bias mea- 
surements by the STAR Beam Energy Scan at RHIC and by the CERES at 
SPS. 



low collision energies, the precise interpretation of data might get difficult due 
to a short plasma phase lifetime, a significant contribution from hadronic phase, 
as well as due to a possible effects of stopping on the bulk evolution [37] . 

Very recently the STAR collaboration completed an even more interesting 
analysis at high energy heavy ion collisions (up to yfs = 200 GeV) following the 
second method mentioned above |22] . By measuring the difference in elliptic flow 
of positive / negative pions for events binned according to their charge asymmetry 
A (at a given centrality class), they found a linear dependence of the splitting 
v 2 ~ v 2 = r e^ 0n the charge asymmetry A. They extracted the slope of this 
linear dependence, that is to be compared with the magnitude of the charge 
quadrupole r e . The comparison between the STAR measured slope parameter 
and our computed charge quadrupole r e shows very good agreement over a 
wide range of centrality for the magnetic field with DH magnitude and lifetime 
of t ~ 5 fm/c or with BS magnitude and lifetime between t = 2.5 — 3 fm/c. 
These measurements at high collision energies and the success of our prediction 
provide a strong evidence for the proposed electric charge quadrupole induced 
by the chiral magnetic wave. 

4 Outlook 

To summarize, we have shown that the chiral magnetic wave induces an elec- 
tric quadrupole deformation of the quark-gluon plasma created in heavy ion 
collisions, and leads to the observable splitting of elliptic flows of positive and 
negative pions. Quantitative estimates of the magnitude and sign of the effect 
(at both low and high collision energies) as well as its linear dependence on 
the charge asymmetry and its centrality dependence (at high collision energies) 
seem to match the experimental data. These results in our opinion can be 
viewed as a strong evidence for the CMW-induced quadrupole deformation of 
the quark-gluon plasma. 

While these developments are very encouraging, more detailed studies still 
have to be done. On the theoretical side, a more sophisticated computation of 
the CMW needs to be performed in an expanding background and taking into 
account the spatial and temporal dependences of the temperature and magnetic 
field. The modeling of the medium and its expansion dynamics can be signif- 
icantly improved by coupling the CMW evolution with the well-developed hy- 
drodynamic simulations. We emphasize that the magnetic field and its lifetime 
are still the main incertitudes that bear significant influence on the quantitative 
results. Progress was made recently on the determination of the electric con- 
ductivity of a pure gluon plasma [29j [30] and we hope that these calculations 
will be extended to the real quark-gluon plasma. Once the electric conductivity 
is known, it would be interesting to perform dedicated simulations of the mag- 
netic field induced inside the plasma. On the experimental side, measurements 
can be developed to further test the proposed V2 splitting. For example, from 
the recently completed U+U collisions, one can look for events with very small 
impact parameter but sizable elliptic flow (the so-called "body-body" collisions) 
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[3T] . These events are believed to have negligible magnetic field and therefore no 
effect from CMW-induced charge quadrupole. Therefore if these events could 
be selected successfully, then one can test the disappearance of the vi splitting. 
One may also test the linear relation of v-x splitting on the charge asymmetry 
A using the data from the Beam Energy Scan. It would also be extremely in- 
teresting to check if this effect, with a linear relation between v-i and A, still 
persists at the LHC energy. 
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